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The doping dependence of the superconducting state for the 2D one-band Hubbard Hamiltonian 
is determined. By using an Eliashberg-type theory, we find that the gap function Ak has a d x 2_ y 2 
symmetry in momentum space and T c becomes maximal for 13 % doping. Since we determine the 
dynamical excitations directly from real frequency axis calculations, we obtain new structures in the 
angular resolved density of states related to the occurrence of shadow states below T c . Explaining 
the anomalous behavior of photoemission and tunneling experiments in the cuprates, we find a 
strong interplay between d-wave superconductivity and dynamical spin fluctuations. 

74.20. Mn,74.72.-h,74.25.Jb 



S3 

33 
in 

On 



> 
N£> 
CO 

NO 

o 
in 

On 



S3 
O 

o 



Despite important progress, the nature of the super- 
conducting pairing mechanismof the High-T c -materials 
is still controversial. Due to their unconventional be- 
havior in the normal as well as in the superconducting 
state various solely electronic pairing mechanism were 
proposed JIJ . During the last years the symmetry of the 
superconducting order parameter in momentum space 
was studied intensively, because it probably holds the key 
for an understanding of the High-T c systems and the role 
of spin fluctuations within the cuprates Recent an- 
gular resolved photoemission (ARPES) experiments and 
phase sensitive measurements of the gap function 
clearly favor a d x 2_ y 2 or an anisotropic s x 2 +y 2 symmetry 
of the gap in momentum space over an isotropic s-wave 
scenario. Moreover the maximum of T c upon doping, the 
occurrence of an additional dip in the ARPES spectra of 
Bi2Sr2CaCu208+<5 found by Dessau et al. || and the dip 
at lo = ±3A in superconductor-insulator-superconductor 
(SIS) tunneling measurements || are important exper- 
iments that may be significant clues for the pairing in- 
teraction. In particular, the recent observation of shad- 
ows of the Fermi surface (FS) in the paramagnetic state 
by Aebi et al. ficj] and their interpretation in terms of 
short-range antiferromagnetic correlations could be the 
Smoking Gun of a spin fluctuation pairing mechanism. 
Therefore, the behavior of the shadow states below T c is 
of great importance for an understanding of the super- 
conductivity in the High-T c systems. 

Theoretically, a favorite model to study a purely elec- 
tronically mediated pairing interaction in the CuC>2 plane 
is the 2D one-band Hubbard Hamiltonian. Recently, 
it was demonstrated within an Eliashberg-type theory 
based on the spin-fluctuation mechanism, that there ex- 
ists a superconducting ground state below T c w 0.02< 
with a d-wave symmetry of the order parameter [[H] ^| . 
Despite these interesting results, the dynamical proper- 
ties and their relation to the strong antiferromagnetic 
correlations are far from being understood, because the 
relevant strong coupling equations were solved on the 



imaginary frequency axis which gives no direct access to 
the dynamical excitation spectrum. A first step to deter- 
mine the excitation spectrum from a real axis calculation 
was achieved in an important study by Dahm et al. |l4| . 

In this Letter, we present results for the excitation 
spectrum in the superconducting state and explain the 
anomalous behavior of the ARPES spectra and tun- 
neling measurements. In addition, we determine the 
doping dependence of the superconducting state and 
investigate the interdependence of the superconduct- 
ing phase and the strong antiferromagnetic correlations 
within the cuprates. Our theory is based on a strong cou- 
pling Eliashberg-type approach for the one-band Hub- 
bard Hamiltonian with nearest neighbor hopping inte- 
gral t = 0.25 eV, bare dispersion £o(k) = —2t[cos(k x ) + 
cos(fcj,)]— fi with chemical potential /i, and local Coulomb 
repulsion U = At. Since we are using our new numerical 
method Esj for the self consistent summation of all bub- 
ble and ladder diagrams (fluctuation exchange approxi- 
mation, FLEX 1 1 6 1 ) on the real frequency axis, we obtain 
directly the quasi particle excitation spectrum below T c . 
The superconducting state is treated in the Nambu for- 
malism where the diagonal and off-diagonal Greens func- 
tion in the matrix notation can be expanded in terms of 
the Pauli matrices fj [i = 0, 1, 3): 



G(k,u) 



ujZ(k,qj)T - (e (k) +x(k,a;))f3 + 0(k,u)fi 
(coZ(k, lo)Y - (e (k) + X (k, - (0(k, W )) 2 



Here, w(l — Z(k, ui)) and %(k, lo) are the diagonal ex- 
pansion coefficients of the electronic self energy matrix, 
whereas </>(k, lo) — A(k, Lo)Z(k, lo) is the anomalous self 
energy and A(k, lo) the gap function JlTj]. Hence, we cal- 
culated these functions by solving the Eliashberg equa- 
tions within the FLEX approximation ]l8| , were the k 
and lo dependent pairing interaction and the Greens func- 
tion G(k, lo) are calculated self consistently Jl9[ . 

In Fig. 1 we show the density of states p(co) in the 
superconducting phase which is always of d x 2_ y 2 pairing 
symmetry for different doping concentrations x = 1 — n 
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and for T — 63 K. Here, n is the occupation number per 
site. For intermediate doping (x = 0.14) a pronounced 
superconducting gap appears. Due to the d-wave pair- 
ing symmetry, states in the gap are clearly visible. For 
larger doping (x — 0.18) the antiferromagnetic correla- 
tions are smaller and consequently the pairing interac- 
tion and the superconductivity becomes weaker, while 
for smaller doping (x = 0.12) the magnetic correlations 
are more dominant and we find a sharp increase of the 
quasi particle scattering. This leads to a loss of metal- 
icity and a weaker superconducting gap. In addition, a 
pseudogap as a precursor of the Mott-Hubbard splitting- 
becomes visible in p(uj). 
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FIG. 1. Density of states p(w) for three doping concentra- 
tions. In the upper inset A(T) is plotted for various doping 
values. The dotted line is the corresponding BCS curve with 
an arbitrary chosen T c — 75 K. In the lower inset the doping 
dependence of T c obtained from A(T C ) = is shown. 

In the upper inset of Fig. 1 we present data for the 
superconducting gap function A(T) estimated from the 
distance between the maxima of p(u>) for different dop- 
ing values. For all doping concentrations the gap opens 
much faster below T c than in the corresponding BCS case 
(dotted lines) and saturates quickly. This is due to the 
shift of spectral weight in the pairing interaction to higher 
energies induced by the opening of the superconducting 
gap and therefore to strong feedback effects [jn]. Fur- 
thermore, we find in agreement with Raman scattering 
measurements that the gap opens more rapidly when x is 
decreased p(]| . Interestingly, the tendency that the gap 
function A(T) saturates at higher values of A(0) when 
the doping is decreased is reversed between x = 0.14 and 
x = 0.13. This behavior reflects the increasing stiffness of 
the system due to the oncoming antiferromagnetic phase 



transition and to the reduction of the charge carrier con- 
centration at the Fermi level resulting from the formation 
of a pseudogap. Therefore, the interplay between antifer- 
romagnetic correlations and superconductivity leads to 
an optimal doping concentration. This is demonstrated 
in the lower inset of Fig. 1, where we present in quali- 
tative agreement with the experimental observation the 
doping dependence of T c becoming maximal for x = 0.13 
and T c = 79 K. 

Recently, we presented corresponding results for the 
normal state jfj| and argued that the observation of 
shadows of the FS shifted by Q = (7r, 7r), by Aebi et 
al. Q is due to strong magnetic correlations that can be 
explained in our framework without a staggered antifer- 
romagnetic moment and with a small correlation length 
of 2.5 lattice spacings. Now, we propose that these 
shadow states also exist below T c and that their intensity 
is even increased compared to the paramagnetic state. 
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FIG. 2. Spectral density p(k, u) for two k values (see 
upper inset), (a) p(k, u) for x = 0.16 and T = 43 K 
(solid line) compared with T — 75 K (dashed line). Note, 
T c (x — 0.16) = 75 K. (b) Solid and dashed lines as in 
(a) and dotted line for x = 0.14 and T = 43 K, where 
T c (x = 0.14) = 78.5 K. The labeled peaks are discussed 
in the text. The lower inset shows the k dependence of the 
scattering rate F(k, uf) for uj = 30 meV where the intensities 
are represented in a linear grey scale. Note, the largest values 
of r(k, u>) are along the main FS and its shadow. 

In Fig. 2 the spectral density p(k,oj) for x — 0.16 and 
T = 43 K is compared with p(k, uj) in the normal state. 
In Fig. 2(a) and 2(b), we plot p(k, uj) for k = (26/32vr, 0) 
close to the FS and for k = (it, 6/327r). For the latter 
momentum one expects in the normal phase the shadow 
states. As was shown in Ref. [|l5| these states are clearly 
visible for smaller doping (x = 0.12). In distinction for 
x = 0.14 only a weak maximum and for x = 0.16 no 
shadow states are observable for this k point above T c . 
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FIG. 3. Quasi particle dispersion below T c along the high 
symmetry line T — * X — * M in the neighborhood of the X 
point. The different symbols represent the dominant contri- 
bution to the peaks in p(k,Lo). 

In Fig. 2(a), the peak which is slightly below the Fermi 
level exhibits a superconducting splitting that is almost 
symmetrical with respect to to = 0. Besides the main 
peaks two dips at lo ss ±40 meV are clearly visible. This 
is related to a strongly enhanced quasi particle decay as 
will be discussed below. 

In Fig. 2(b) the main peak in the normal state is 
well above the Fermi level. However, a superconduct- 
ing splitting can be seen below T c . Pronounced quasi 
particle decay processes in analogy to Fig. 2(a) and 
the superconducting splitting of the main band tail for 
\io\ < |Afc| lead to the significant dip structure for 
co > 0. Most interestingly, these structures can only be 
found for k states near the FS and the corresponding 
shadow states, if A& ^ 0. This results from an anoma- 
lous k dependence of the quasi particle scattering rate 
r(k, lo) — Im(ojZ(k, w))/ReZ (k, lo) shown in the lower 
inset of Fig. 2 for lo — 30 meV. Its maxima occur close to 
the FS but also at the position of its shadow with almost 
equal intensities leading to a suppression of the spectral 
weight. We expect this result to be also of importance 
for the unusual transport properties of the cuprates. In 
addition, it is highly interesting that the dip structure in 
Fig. 2(b) is in good agreement with the corresponding 
ARPES measurements by Dessau et al. in the BISCO 
system P,pT|- By decreasing the doping to x = 0.14 the 
distance between the quasi particle peak in the normal 
state and the Fermi level increases. Therefore, the spec- 
tral weight of the superconducting quasi particles in peak 
II decreases compared to a; = 0.16, which one also would 



expect for the superconducting image peak /. However, 
one clearly sees an increase of peak /. In order to show 
that this observation is due to the formation of super- 
conducting quasi particles within the shadow states, we 
analyzed the quasi particle dispersion in the supercon- 
ducting phase. 

In Fig. 3, we present this dispersion e(k) obtained 
from the maxima of p(k, lo). The filled squares indicate 
the main-band whereas the open squares are due to its 
superconducting splitting. Note that both exhibit a flat 
dispersion within a large k range. This is due to the flat 
bands in the paramagnetic state and the fact that only 
the spectral weight below \lo\ « |Ak| is splitted in the su- 
perconducting phase. Furthermore, the states indicated 
by open triangles are the shadow states shifted by Q with 
respect to the main band. In the paramagnetic state and 
for x — 0.16 these shadow bands are also existent, but 
their intensities are much weaker. Therefore, a clear am- 
plification of the shadows below T c occurs. Moreover, we 
indicate with crosses the maxima caused by the formation 
of the dip structure. They result, similar to the shadow 
states, from strong quasi particle decay. For that reason 
it is tempting to assume that these structures are the 
continuation of the shadow band structure below T c for 
a system with purely dynamical antiferromagnetic cor- 
relations. The four band splitting at the FS is due to 
the fact that the frequency dependence of the quasi par- 
ticle decay rate of a given band determines that of its 
superconducting image state. Since both peaks behave 
symmetrically, we find the significant four band splitting 
at the FS. Finally, this amplification of the shadow states 
explains the anomalous doping dependence of the spec- 
tral density shown in Fig. 2(b). Although the distance of 
the peak II and the Fermi level increases upon doping, 
the crossing point of the shadow states and the Fermi 
level is shifted towards k = (it, 6/327r). The anomalous 
behavior of peak / results from an increase of its shadow 
state contribution for smaller doping. 

From the effective interaction of the quasi particles, 
we find that the antiferromagnetic coupling in the super- 
conducting state is even stronger compared to the situa- 
tion above T c . Although the magnetic correlation length 
(k, 3 lattice spacings) increases by less that one lattice 
spacings compared to the paramagnetic phase, a clear 
increase of the effective interaction (« 15%) occurs. Due 
to the retarded coupling of singlet Cooper pairs between 
nearest neighbor sites, the antiferromagnetic correlations 
are stabilized on the corresponding time scale. These re- 
sults indicate not only that the d-wave superconductivity 
is caused by antiferromagnetic correlations, but that the 
presence of the superconducting state is also enhancing 
the short range antiferromagnetic fluctuations. In order 
to demonstrate further consequences of our theory we 
also calculated the tunneling conductance for SIN and 
SIS junctions, because SIS in-plane tunneling measure- 
ments exhibit an interesting dip structure at lo = ±3A 
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whose temperature dependence is much stronger com- 
pared to the main gap. 
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FIG. 4. SIS tunneling conductance a s (eV) for x = 0.16 
normalized by the paramagnetic a n (eV) for different tem- 
peratures. Note the strong temperature dependence of the 
additional dip structure at ±3A. 

In Fig. 4 we present results for the SIS in-plane tun- 
neling conductance a s (eV) — dI(eV)/dV |j| using 

/oo 
dujp(uj)p{LU + eV)(f{uj) - f(u + eV)) 
-oo 

for the tunneling current. Here, e is the unity charge, 
f(uj) the Fermi function, and V the applied bias volt- 
age. By inserting our self-consistently calculated p(uj), 
we find an excellent agreement for both the temperature 
dependence and the position of the dip with the experi- 
mental data. Note that in agreement with SIN measure- 
ments, no pronounced dip in the corresponding tunneling 
can be observed (not shown). It is important to point 
out that our interpretation of the 3A dip structure in 
terms of a strong increase of the scattering rate r(k, uj) 
at the shadow states is in contrast to the explanation 
of Ref. |^2| . Therein its authors argued that the dip in 
the SIS measurements is caused by a sudden increase of 
r(k, uj) at wave vectors k x = ik y where the gap vanishes. 

In conclusion, we calculated the superconducting prop- 
erties of the 2D Hubbard model. We found that the 
order parameter has a d x 2_ y 2 symmetry in momentum 
space and that T c is maximal for the doping concentra- 
tion x = 0.13. It was demonstrated for the first time that 
the shadow states observed by Aebi [|l0| are still present 
and even enhanced below T c . We obtain a four-fold split- 
ting of the quasi particle dispersion at the FS and a 
dominant dip structure in the spectral density p(k, uj) 
next to the shadows of the FS. These new structures are 
caused by purely dynamical antiferromagnetic spin fluc- 
tuations. In view of the recent phenomenological model 
for a spin fluctuation induced anisotropic s-wave pairing 
symmetry in bilayers |23j| we believe that in our theory a 



similar interplay of antiferromagnetic and superconduct- 
ing excitations occurs for multi-layer systems too. From 
our spectral density we calculate the in-plane SIS tun- 
neling conductance and demonstrate that the dip in the 
ARPES and the 3A dip in the tunneling characteristics 
arise for the same physical reason. Therefore, we believe 
that these important experiments and the observation 
of shadows of the FS by Aebi et al. can be explained 
within our theory. All this sheds new light on the super- 
conducting pairing mechanism and supports the role of 
spin fluctuations for the high T c in the cuprates. 
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